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ABSTRACT 

Most stellar remnants so far have been found in binary systems, where they interact 
with matter from their companions. Isolated neutron stars and black holes are difficult 
to find as they are dark, yet they are predicted to exist in our Galaxy in vast numbers. 

We explored the OGLE-III database of 150 million objects observed in years 
2001-2009 and found 59 microlensing events exhibiting a parallax effect due to the 
Earth’s motion around the Sun. Combining parallax and brightness measurements 
from microlensing light curves with expected proper motions in the Milky Way, we 
identified 13 microlensing events which are consistent with having a white dwarf, 
neutron star or a black hole lens and we estimated their masses and distances. The 
most massive of our black hole candidates has 9.3 Mq and is at a distance of 2.4 kpc. 
The distribution of masses of our candidates indicates a continuum in mass distribution 
with no mass gap between neutron stars and black holes. We also present predictions 
on how such events will be observed by the astrometric Gaia mission. 
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1 INTRODUCTION 

Dark stellar remnants, namely neutron stars (NS) and black 
holes (BH), are difficult to study as they are generally very 
hard to find. Yet, they are important ingredients in our un¬ 
derstanding of stellar evolution, mass distribution, galaxy 
evolution and structure, dark matter balance, etc. 

Both NS and BH are typically being discovered in bi¬ 
nary systems thanks to interactions with a companion star 
(e.g., Ziolkowski 2010, Bachetti et al. 2014). Single neutron 
stars can also be found as radio pulsars when their align¬ 
ment is fortunate enough so that their radio beams point 
toward us (e.g.,Posselt et al. 2008, Popov 2012), as well as 
in gamma rays (e.g., Abdo et al. 2010). Masses of NS and 
BH have been measured only in a small few dozen cases and 
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their mass functions has been estimated (e.g., Kiziltan et al. 
2013, Ozel et al. 2010, Kochanek 2015). The observational 
data collected so far indicate a clear gap between neutron 
stars (with masses up to 2 M©) and black holes (masses 
from 6 Mq), e.g.,Bailyn et al. 1998. Interestingly, the the¬ 
oretical predictions of the end products of stellar evolution 
can not reproduce the gap and rather suggest a continuum 
in mass distribution of remnants (Fryer & Kalogera 2001). 
This puzzle remains unresolved and the larger new sample 
of isolated NS and BH could unveil the answer about the 
mass gap. 

Gravitational microlensing is the only method capa¬ 
ble of finding isolated stellar-mass remnants and deriving 
their masses and the mass function (e.g., Paczynski 1986, 
Dai et al. 2015), as it does not require light emitted by the 
lens to detect it (Paczynski 2003). Gould (2000b) estimated 
that about 4 per cent of lenses towards the Galactic Bulge 
should be due to dark remnants: neutron stars (~3 per cent) 
and black holes (~ 1 per cent). There should also be about 17 
per cent of luminous, yet still very faint, white dwarf lenses. 
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More recently, Oslowski et al. (2008) simulated the BH and 
NS microlensing events in the Milky Way and estimated 
that current microlensing surveys should expect from 4 to 
10 such events every year, depending on the model. Cur¬ 
rently, OGLE-IV (Udalski et al. 2015) detects more than 
2000 microlensing events every year, therefore there should 
be at least a dozen of lensing black holes among them, yet 
none were conclusively identified as such. The primary rea¬ 
son for this is that the vast majority of microlensing events 
are standard single point lens events (e.g., Paczynski 1996, 
Wyrzykowski et al. 2015) and such events do not allow for a 
determination of lens mass and distance. Additional effects 
and observations are required in order to break the degen¬ 
eracies and to compute the mass of the lens. According to 
Gould (2000a), to measure the mass of the lens only two 
additional observables are needed, 6e and tte: 

M=— (1) 

K-TTe 

where Oe is the angular Einstein radius of the lens, con¬ 
stant K = 4:G/{c^ AU) = 8.144 mas/MQ and tte is microlens¬ 
ing parallax. The measurement of 6e is possible in very rare 
cases, when the source star disk is being resolved during 
a high magnification or caustic crossing event and can be 
used as an angular ruler (finite source effect) (e.g., Zub et al. 
2011, Lee et al. 2010). Another way is to detect both source 
and lens and to measure the relative proper motion, which 
leads to 0 e (e.g., Gould et al. 2004, Kozlowski et al. 2007, 
Bennett et al. 2015, Batista et al. 2015). Yet another pos¬ 
sibility is via precise astrometric measurements of the cen¬ 
troid motion during an event, however, this has not yet been 
achieved as the displacements to be measured are of order 
of 1 mas and are challenging for current facilities. In the 
near future, however, such measurements might be possible 
with the Gaia mission (e.g., Belokurov & Evans 2002, Proft 
et al. 2011, Wyrzykowski & Hodgkin 2012) and WFIRST 
(The Wide-Field Infrared Survey Telescope, Bennett & Rhie 
2002, Spergel et al. 2015,Gould & Yee 2014). 

Obtaining the second parameter from Eq. 1, tte, is 
somewhat easier. Microlensing parallax is an analog of stan¬ 
dard parallax, and relies on detecting differences in an mi¬ 
crolensing event caused by different viewing angles. Re¬ 
cently, a very successful programme of observing selected mi¬ 
crolensing events from the Spitzer satellite, separated from 
the Earth more than 1 AU, has yielded parallax measure¬ 
ments for otherwise standard events (Yee et al. 2015). 

However, in the case of long timescale events, the Earth 
changes its position significantly in its orbit around the Sun 
and therefore Earth-based parallax can be measured using 
only ground-based observations (Smith et al. 2005, Gould 
2004). Long time scale events are also more likely to be 
owing to high mass lenses. In an event showing the effect of 
parallax, the observers motion changes the overall shape of 
the microlensing light curve, which deviates from a standard 
Paczynski curve. 

Measuring tte, combined with certain assumptions, 
e.g., on the proper motions of the lens and the source, pro¬ 
vides probability distributions for the masses and distances 
of lenses. The terrestrial parallax effect has been used to 
suggest that three events from the MACHO survey are can¬ 
didates for isolated black-hole lenses (Bennett et al. 2002), 
however, none of them was ever confirmed, despite various 


attempts including X-ray follow-up (e.g., Maeda et al. 2005, 
Nucita et al. 2006). Mao et al. (2002) also used the par¬ 
allax effect to constrain the mass of the candidate black- 
hole lens to about 10 Mq in an event from the OGLE- 
H survey (OGLE-1999-BUL-32), which was also the event 
with the longest Einstein radius crossing time found to date 
(tB=640 days). However, X-ray follow-up only placed an up¬ 
per limit on the X-ray luminosity (Revnivtsev & Sunyaev 
2002). There was also a candidate for a binary black hole 
system, OGLE-2005-SMC-001, where parallax and weak bi¬ 
nary signatures were used to conclude that the components 
had masses 7 Mq and 3 Mq (Dong et al. 2007). 

In this paper we present a comprehensive study of the 
ground-based parallax microlensing events found in the data 
of the OGLE-HI project (section 2). We search for new can¬ 
didates for dark remnant lenses (section 3), including can¬ 
didates for black-holes and derive the mass function for iso¬ 
lated neutron stars and black-holes (section 4). We discuss 
the results in section 5 and conclude in section 6. 


2 DATA 

The photometric data used here are the same as those 
used for our search for standard microlensing events 
(Wyrzykowski et al. 2015). In brief, the Optical Gravita¬ 
tional Lensing Experiment in its third generation (OGLE- 
HI, Udalski et al. 2008) operated from 2001 until 2009 and 
used the 1.3m Warsaw Telescope, located at the Las Cam- 
panas Observatory, Chile, operated by the Carnegie Institu¬ 
tion for Science. OGLE-HI observed the Galactic Bulge in 
177 fields; 0.34 sq. deg each over 8 mosaic CCDs. We selected 
the best observed central 91 fields covering, in total, 31 sq. 
deg and containing about 150 million sources, each having at 
least 250 observations. The survey was conducted primarily 
in the /-band, calibrated to the Cousins /-band (Szymahski 
et al. 2011), and in the dense Bulge regions reached down 
to about 20.5 mag. The typical sampling was, on average, 
once per three nights, with the exception of the most central 
fields, which (from 2005) were observed about three times 
per night, yielding up to 2500 data points for those regions. 
Fig. 1 shows a map of the OGLE-HI fields with all microlens¬ 
ing events exhibiting a parallax effect marked. 

U-band images were collected occasionally, with up to 
35 data points available per light curve, providing an oppor¬ 
tunity to derive average colours of the stars. 

The images in both bands were reduced on-the-fly soon 
after the observations using Difference Imaging Analysis 
(DIA, Wozniak 2000), which allowed for real-time discover¬ 
ies of on-going microlensing events via OGLE’s Early Warn¬ 
ing System (EWS\ Udalski 2003). At the end of OGLE-HI 
the data were re-reduced using better quality reference im¬ 
ages (Udalski et al. 2008). This is the primary database in 
which we conducted our search. The microlensing events dis¬ 
covered among the 150 million light curves in the database 
were then additionally re-reduced further, using more pre¬ 
cise light centroid positions, in order to improve their pho¬ 
tometry. All microlensing modelling was performed using 


^ http://ogle.astrouw.edu.pl/ogle3/ews/ews.html 
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Figure 1. Map of the OGLE-III bulge fields selected for search for parallax events (red squares). Black squares show other OGLE-III 
fields with not enough observations. Circles show 59 parallax microlensing events, while filled circles mark candidates for dark remnants 
lens. The background image was taken by Krzysztof Ulaczyk. 


this new photometry, but the search itself and the simula¬ 
tion of events used the main data base, in order to main¬ 
tain consistency. The photometric error-bars were corrected 
using scaling coefficients derived for each field/chip (see 
Wyrzykowski et af. 2009). A few frames were removed from 
the anafysis as they were taken under very poor conditions 
(mostly from the very first season of OGLE-ffI, when the 
survey was just starting) and obvious outlying data points 
were removed using a sigma-clipping method. 

3 SEARCH FOR REMNANT LENSES 

We explored the vast OGLE-fff database of light curves con¬ 
taining about 150 million sources over a number of stages in 
order to find microlensing events showing a parallax effect, 
which coufd be attributed to stellar remnant lenses. 

3.1 Preselection of objects 

As the first step we take the simple and quick to compute von 
Neumann (rj) and skewness ( 7 ) statistics, previously com¬ 
puted for all stars in the 91 OGLE-fff fiefds in the search 
for dips in the OGLE fight curves (Rattenbury et al. 2015). 
The combination of those two statistics is sensitive to tem¬ 
poral changes in brightness - a negative value of skewness 
indicates a brightening, like in microlensing events. We then 
visually selected a couple of dozen long-lasting microlens¬ 
ing events (found during the search for standard events 
in Wyrzykowski et al. 2015) and used them to derive the 
boundary in the rj and 7 parameter space. 

The cut was able to limit a large sample of stars to 
only those which exhibited a long-term smooth variations 
(16,060). Among those we removed stars which were identi¬ 
fied as long period variable stars in Soszyriski et al. (2013) 


and ghost events Wyrzykowski et al. (2015), which often 
were multiplied occurrences of the same object. We were 
then left with just 3845 light curves. Those were then fitted 
with the standard Paczynski microlensing model for a sin¬ 
gle lens, with a blending parameter fixed at a set of values 
ranging from 1 to 0 . 1 , to avoid models not converging numer¬ 
ically. At this stage we used a simple and quick to compute 
Levenberg-Marquardt Least-Square (LSQ) minimization al¬ 
gorithm in order to filter out genuine microlensing events 
with parallax from other outburst-like and variable objects. 
The parameters of the standard model were the following: 
to* - time of the closest lens-source approach; u'q - minimum 
impact parameter; - Einstein radius crossing time (the 
events time-scale); /q* - baseline magnitude; /|* - blending 
parameter (defined as the ratio of the source flux contribu¬ 
tion to the total baseline flux). 

In the next step, the parameters of the standard model 
were used to seed the search for a parallax microlensing 
model (Gould 2004). The parameters of the parallax model, 
computed in the geocentric frame, included five parameters 
as in the standard model (but indexed with tt in order to dis¬ 
tinguish from their standard model equivalents) and two ad¬ 
ditional parameters: tten and ttee, the North and East com¬ 
ponents, respectively, of the parallax vector, tte. The geo¬ 
centric parallax model also requires a reference time (tg***)! 
which was chosen to be tg*- 

While fitting the parallax model with LSQ the blending 
parameter (/s) was again fixed at set values from 1 to 0.1 in 
order to avoid non-convergence of models (but independent 
on fs derived for the standard model, as the parallax effect 
can sometimes mimic blending in the standard model). 

We then computed the difference between the good¬ 
ness of fit between the parallax and standard models (Ax^), 
which was used as the main indicator for the significant par- 
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allax signal. We also added an empirical parameter, ampl, 
which was the amplitude of the light curve, in order to 
exclude very shallow and long-term variable objects, not 
caused by microlensing. 

We have visually inspected all 3845 objects remain¬ 
ing after the rejection of variable stars and ghosts and se¬ 
lected 59 microlensing events which were clearly exhibiting 
parallax-like deviation from the standard model. The vi¬ 
sually selected sample (hereafter referred to as the golden 
sample) allowed us to construct a set of preliminary cuts 
on various parameters of the fitted models in order to nar¬ 
row down and automatise the selection. The series of cuts 
used microlensing parameters from both standard and par¬ 
allax LSQ models, along with quantities like the amplitude 
of variation and other statistics. The entire sequence of the 
cuts applied is described in Table 1. The results of the cuts 
is a set of 217 objects, which includes all 59 golden sample 
parallax events. 

3.1.1 Random Forest classifier 

The cuts on parameters described above were not ideal as 
they were simple linear cuts through the parameter space 
and the selected objects still required a visual inspection. In 
order to prepare a fully automated search pipeline, it was 
necessary to replace the visual inspection with more com¬ 
plicated selection criteria, which can be found using Ma¬ 
chine Learning algorithms. We therefore substituted the hu¬ 
man with a similarly skilled Random Forest (RF) classiher 
(Breiman 2001, Wyrzykowski et al. 2015). 

The training of the RF was performed on the entire sam¬ 
ple of 217 events which remained after the cuts, using their 
parameters, e.g., microlensing fit parameters of standard and 
parallax model, rj, 7 , Ax^ between the goodness of fit be¬ 
tween parallax and standard models. The cross-validation of 
that training (as run in the Weka^ package) yielded about 
81 per cent correct classihcations, however, the entire golden 
sample was retrieved without any losses when run on the 
whole sample. This guarantees that the classifier did not 
exclude good events, but performed well on removing bad 
cases, i.e.,with high completeness. The classifier would be 
used for simulated events only, hence its primary role was 
not to detect and remove non-microlensing events, but to 
recognise microlensing events with the parallax effect. 

Table 2 lists all 59 golden sample parallax microlensing 
events, their names, OGLE database IDs, coordinates and 
length, dehned as the time from detectable microlensing de¬ 
viation computed from the parallax model (i.e.,by Icr from 
the baseline). 

3.1.2 MCMC Modelling 

Each of the selected events from the golden sample was mod¬ 
elled with the Markov chain Monte Carlo (MCMC) method. 
Full MCMC runs were made across on an extended range of 
parameters allowing all degenerate solutions for microlens¬ 
ing parallax model to be found (e.g., Gould 2004, Smith et al. 
2005) as well as obtaining probability distributions for each 
of the parameters of the model. To find the microlensing 

^ http;//www.cs.waikato.ac.nz/ml/weka/ 


parameters and their 1-a limits we use the Affine Invariant 
MCMC Ensemble sampler - emcee (Foreman-Mackey et al. 
(2013)). Note, in the MCMC parallax microlensing models 
we allowed for a small amount of negative “third light”, so 
called, negative blending. In our definition of blending, the 
negative blending occurrs when /s > 1 (Smith et al. 2007). 
Small amount of negative blending is an indication of possi¬ 
ble over-estimation of the background level on the reference 
image, locally influenced by numerous unresolved stars (see 
also Wyrzykowski et al. 2015). We also tested the conver¬ 
gence of the MCMC chains for unconstrained /s and found 
that some of the solutions were nonphysical (/s > 1 ), hence 
those were removed. 

We denote the microlensing parallax model parameters 
obtained from the MCMC modelling with no additional in¬ 
dices, in order to distinguish them from the parameters from 
the quick LSQ modelling performed earlier in the search 
pipeline. 

The parameters from the MCMC modelling for all 59 
events are available on-line® and in the on-line version of 
the paper. Table 3 shows the microlensing parallax model 
parameters for 13 events selected as the most probable can¬ 
didates for having dark stellar remnants as lenses (see next 
section). 


3.1.3 Colour-magnitude diagram selection 

In order to limit our sample of events to those with sources 
at around 8 kpc (e.g., Pietrukowicz et al. 2014), i.e.,in the 
Galactic Bulge, we constructed a colour-magnitude diagram 
(CMD), as shown in Fig. 2. The 7-band magnitudes of the 
sources were de-blended using Jq and fs from the microlens¬ 
ing parallax MCMC models and were corrected for extinc¬ 
tion, based on Nataf et al. (2013) extinction maps for OGLE- 
III. We must note here, that due to the scarcity of R-band 
observations during the OGLE-III phase, the colour infor¬ 
mation was, in the case of most events, only available for 
the baseline. Hence our V — I colour is only an approx¬ 
imate colour of the source. However, as already noted in 
Wyrzykowski et al. (2015), the bulge Red Clump Giants 
(RCGs) are typically brighter than most of other stars on 
the CMD, hence even if blended with a Main Sequence star 
in the baseline, a RCG source would not change position 
greatly on the CMD. 

We defined the Red Clump Region (RCR) as 0.5 < 
{V — 7)0 < 2.5 and 12 < ( 7 source)o < 16 for extinction-free 
magnitudes - this selected 26 sources. For sources from the 
RCR we assumed their distance to be Ds = 8 kpc. 

The source in one of our events, PAR-42, has been 
observed spectroscopically with the 2.5m Irenee du Pont 
telescope at Las Campanas Observatory in April 2014 (see 
Pietrukowicz et al. 2015 for details). The low-resolution 
spectra, obtained with the B&C spectrograph, clearly 
showed that the source is a K3 Red Clump Giant, confirming 
our conclusions based solely on the CMD from OGLE-III. 


^ http://ogle.astrouw.edu.pl 



BH, NS and WD Candidates from Microlensing 5 


Table 1. Cuts applied to the entire OGLE-III database in order to detect events with a clear parallax signal. 


Cut number 

Formula 

Descriptio: 




N objects remaining 

0 

91 fields, at least 250 datapoints 

initial database 



150X 10® 

1. 

< 0 log(—> -1.25{log(l/T7) + 0.75 

up-wards £ 

ind long-term deviations 



16060 

2. 

NOT var &c NOT ghost 

removed known variable stars and ghost ev 

ents 


3845 

3. 

2050 < JD(tQ) - 2450000 < 5000 

maximum 

within the data range 



3645 

4. 

> 50 & Ax^ > 75(xJ/Ndof - 3) S- 60 

parallax ip 

Lodel significantly better than st£ 

indard 


679 

5. 

fg > 20 & f|? < 650 <S4 tg > 13{xJ/Ndof - 3) + 5 

non-parall. 

ax time-scale (duration of the ev 

ent in days) abo'v 

•e thresholds 

389 

6. 

ampl > 0.16 mag 

significant 

amplitude of variation 



318 

7. 

> 50 d 

parallax m 

Lodel time-scale long enough 



245 

8. 

xJ/Ndof < 10 

parallax m 

Lodel well fit 



217 

(9.) 


visual inspection/Random Forest classifier 



59 


12 



^”-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 

{V - /)o [mag] 


Figure 2. Extinction free colour-magnitude diagram for sources 
in our sample of parallax events from OGLE-III. The back¬ 
ground is the extinction-corrected CMD of stars in the sub-field 
BLG205.3. The green box indicates our Red Glump Region. Dark 
remnant events with probability above 75 per cent are marked in 
red. 

3.2 Simulations and Detection Efficiency 

Our experiment is naturally not sensitive to every parallax 
event. For example we might have missed an anomaly due to 
parallax in an event during a gap in our observation season. 
In order to correct for the detection efficiency we conducted 
simulations of parallax microlensing events. 

The events were simulated in a similar fashion as de¬ 
scribed in Wyrzykowski et al. (2015). Mock microlensing 
events were generated on top of light curves from the OGLE- 
III Bulge data base, preserving any intrinsic variability and 
original noise properties originating from varying observ¬ 
ing conditions (seeing, airmass, background). Thanks to 
drawing random stars from the OGLE data base we re¬ 
produced the underlying observed luminosity function for 
the baseline of the simulated events. Blending parameters 
(/s) were drawn from an empirical distribution obtained 
based on the comparison between observed OGLE stellar 
density and archival Hubble Space Telescope (HST) /-band 
images (as shown in Wyrzykowski et al. 2015). The im¬ 
pact parameter (uq) was taken from a flat distribution from 
0 to 1 and the time of minimum approach (to) was ran¬ 



Figure 3. Detection efficiency for parallax events as a function 
of and TTE (left panel). Dots denote detected microlensing 

parallax events. Right panel shows the efficiency as a function of 
a mass and distance to the lens, for assumed, fixed proper motion 
of = 5 mas/yr. 

domly chosen from the entire range of the covered epochs 
(2100 < JD — 2450000 < 5000). We simulated events with 
time-scales Ie ranging from 1 to 1000 days, and parallaxes 
tte ranging from 0.0001 to 0.5, with a random angle for 
tte, decomposed into North and East components (ttbn and 
ttee). 

After applying all cuts as for the real data (including 
Random Forest classifier), the efficiency was derived as a 
function of two parameters: Ie and tte. 

Figure 3 shows the derived detection efficiency as a func¬ 
tion of Ie and tte. All parallax microlensing events detected 
in the OGLE-III data are marked as dots. When combining 
both parameters, the efficiency can be also expressed as a 
function of plausible mass of the lens (right panel in Figure 
3), when assuming a fixed value of prei = 5 mas/yr, the most 
likely value of relative proper motion (see next section). 

As expected, we are not likely to detect the parallax 
signal in short time-scale events with small parallax mod¬ 
ulation, hence, for low-mass lenses or lenses located in the 
Galactic bulge. Decreasing efficiency for very long events 
and with very strong parallax is due to the fact that such 
events usually exceed the data span we covered in OGLE- 
III and their parallax modulations were so strong that there 
was hardly any baseline in the light curve to conclude on its 
microlensing nature. Our search for parallax events has the 
highest efficiency for long events or large parallaxes, reaching 
its maximum of about 55 per cent for time-scales of about 
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Table 2. Parallax microlensing events found in OGLE-III. Columns show the unique name of the lens, OGLE database id, right-ascension, 
declination, OGLE-III EWS match and length of the light curve. 


Lens name 

OGLE3-ULENS- 

OGLE DB ID 
[field, chip, star no] 

R-Aj2000 

[h:m:s] 

DECj2000 

[d:m:s] 

EWS 

OGLE- 

length* 

[days] 

PAR-01 

BLG196.1.118200 

18:04:48.79 

-29:40:31.1 

_ 

2840.8 

PAR-02 

BLG205.3.159237 

17:57:23.14 

-28:46:32.0 

2006-BLG-095 

2306.4 

PAR-03 

BLG242.7.31173 

18:08:0.42 

-26:39:15.4 

— 

2240.3 

PAR-04 

BLG342.5.73806 

17:46:25.67 

-22:51:37.3 

— 

2225.4 

PAR-05 

BLG138.1.192949 

17:46:35.45 

-33:46:19.8 

2004-BLG-361 

1947.5 

PAR-06 

BLG225.8.163880 

18:03:52.40 

-27:54:1.3 

2006-BLG-393 

1566.3 

PAR-07 

BLG102.7.44461 

17:55:59.44 

-29:38:28.8 

2005-BLG-474 

1244.8 

PAR-08 

BLG130.1.122103 

17:47:34.77 

-34:25:50.3 

2005-BLG-351 

1223.2 

PAR-09 

BLG172.8.614 

17:54:11.82 

-31:31:57.9 

— 

1141.5 

PAR-10 

BLG156.8.6S730 

17:54:4.49 

-32:40:38.2 

— 

1044.2 

PAR-11 

BLG197.3.171180 

18:07:33.38 

-29:20:16.2 

2003-BLG-263 

1030.3 

PAR-12 

BLG129.6.16369 

17:42:58.78 

-34:11:13.2 

2005-BLG-059 

991.0 

PAR-13 

BLG333.2.53483 

17:35:55.95 

-27:16:2.1 

2002-BLG-061 

972.9 

PAR-14 

BLGIOI.7.52163 

17:53:28.47 

-29:57:44.7 

2008-BLG-318 

970.6 

PAR-15 

BLG223.5.115009 

17:58:27.13 

-27:26:0.9 

2008-BLG-096 

935.2 

PAR-16 

BLG172.5.160833 

17:54:55.05 

-31:00:40.5 

2003-BLG-032 

928.9 

PAR-17 

BLG249.2.92350 

18:06:6.30 

-25:59:12.8 

2008-BLG-157 

927.8 

PAR-18 

BLG205.5.12205 

17:56:7.08 

-28:41:37.3 

— 

848.8 

PAR-19 

BLG194.3.95955 

17:51:48.90 

-29:17:52.9 

2005-BLG-372 

826.8 

PAR-20 

BLG197.6.9S763 

18:06:59.34 

-29:23:39.7 

2008-BLG-098 

800.9 

PAR-21 

BLG122.1.184151 

17:50:4.59 

-34:58:15.3 

2006-BLG-005 

790.3 

PAR-22 

BLG158.5.38090 

18:00:20.30 

-32:15:11.3 

2002-BLG-334 

755.1 

PAR-23 

BLG224.3.40806 

18:02:53.14 

-27:41:15.7 

— 

747.5 

PAR-24 

BLG236.7.59423 

18:09:13.92 

-27:12:12.6 

— 

747.2 

PAR-25 

BLG251.4.149142 

18:11:29.35 

-25:38:24.0 

2005-BLG-036 

747.1 

PAR-26 

BLG130.4.125656 

17:47:27.03 

-33:58:8.4 

2006-BLG-020 

745.9 

PAR-27 

BLG234.6.218982 

18:04:45.71 

-26:59:15.3 

2005-BLG-086 

735.4 

PAR-28 

BLG197.1.69448 

18:07:41.45 

-29:42:39.1 

2005-BLG-020 

716.7 

PAR-29 

BLG218.3.139976 

18:08:46.04 

-28:12:34.6 

2006-BLG-031 

680.6 

PAR-30 

BLG104.1.148712 

17:59:44.30 

-29:41:7.3 

2008-BLG-545 

677.8 

PAR-31 

BLG167.8.67944 

18:02:41.84 

-32:04:39.6 

— 

672.0 

PAR-32 

BLG195.2.8006 

17:53:55.65 

-29:22:41.1 

2008-BLG-223 

667.1 

PAR-33 

BLG208.5.131806 

18:04:4.02 

-28:38:38.8 

2006-BLG-251 

661.0 

PAR-34 

BLG180.6.142446 

17:51:38.96 

-30:32:16.9 

— 

604.5 

PAR-35 

BLG354.3.45010 

17:35:25.94 

-23:31:0.7 

2003-BLG-297 

583.9 

PAR-36 

BLG208.3.222797 

18:06:21.76 

-28:44:46.3 

— 

555.0 

PAR-37 

BLG195.5.141274 

17:53:16.57 

-28:58:49.8 

— 

552.4 

PAR-38 

BLG155.8.123452 

17:51:25.67 

-32:38:6.1 

— 

510.9 

PAR-39 

BLG171.4.87121 

17:53:34.61 

-31:01:12.3 

2008-BLG-013 

505.2 

PAR-40 

BLG157.3.1531 

17:57:54.73 

-32:25:31.9 

2005-BLG-454 

492.5 

PAR-41 

BLG226.1.1251 

18:07:8.81 

-27:59:22.5 

2003-BLG-459 

467.3 

PAR-42 

BLG134.5.193547 

17:56:59.79 

-33:55:15.6 

2005-BLG-061 

463.1 

PAR-43 

BLG182.3.104094 

17:58:30.91 

-30:34:7.6 

2009-BLG-049 

459.7 

PAR-44 

BLG180.7.173170 

17:51:21.97 

-30:41:1.2 

2007-BLG-282/2007-BLG-303 

452.8 

PAR-45 

BLG226.8.198790 

18:05:51.80 

-27:52:23.3 

2006-BLG-015 

428.5 

PAR-46 

BLG156.2.110020 

17:55:5.17 

-32:28:51.4 

2006-BLG-023 

427.9 

PAR-47 

BLG179.2.183955 

17:51:21.96 

-30:41:1.2 

2007-BLG-282/2007-BLG-303 

417.3 

PAR-48 

BLG197.1.101151 

18:07:8.11 

-29:39:34.0 

2008-BLG-375 

410.8 

PAR-49 

BLG185.7.91655 

18:05:34.74 

-30:44:20.2 

2007-BLG-035 

398.1 

PAR-50 

BLG227.7.16177 

18:08:41.73 

-27:49:14.9 

— 

338.3 

PAR-51 

BLG249.2.60315 

18:06:34.67 

-26:01:16.0 

2003-BLG-175 

333.6 

PAR-52 

BLG241.2.105532 

18:07:5.68 

-26:36:33.1 

2008-BLG-038 

325.6 

PAR-53 

BLG195.2.18171 

17:54:5.53 

-29:23:11.0 

2004-BLG-070 

323.3 

PAR-54 

BLG104.7.157693 

17:58:29.56 

-29:30:54.3 

2006-BLG-366 

322.4 

PAR-55 

BLG138.7.140195 

17:45:11.29 

-33:38:55.2 

2006-BLG-061 

318.2 

PAR-56 

BLG155.1.109582 

17:52:23.34 

-32:37:57.9 

— 

296.5 

PAR-57 

BLG241.6.137439 

18:05:24.43 

-26:25:18.8 

2007-BLG-349 

261.5 

PAR-58 

BLG188.5.139006 

17:57:58.60 

-29:48:48.4 

— 

216.7 

PAR-59 

BLG194.7.63636 

17:49:44.34 

-29:28:34.3 

— 

115.4 


* Length is the duration of the entire event, defined as the time between measurable microlensing deviation, depending on the photometric error-bar (see 

text). 


300 days and parallaxes ~ 0.1. There is also an extension 
at very long time-scales (~ 600 days) for smaller tte ~ 0 .1. 
In terms of approximate masses and distances, our pipeline 
is most sensitive for lenses heavier than about 0.7 M© and 
at a distance between 1 and 4 kpc. Our detection efficiency 
also remains high (more than 50 per cent) for very heavy 
lenses of about 10 Mq, where the simulation for efficiency 
stopped. 


4 DARK REMNANT CANDIDATES 


For a regular standard microlensing event it is typically im¬ 
possible to obtain a unique solution for the distance and 
mass of the lens, since there is only one physical parameter 
(tfi) obtained from the light curve model, which is linked 
with source and lens distance, lens mass and lens-source 
relative velocity. However, for parallax events the situation 
improves as we additionally measure the microlensing par- 
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Table 3. Microlensing parallax model parameters as obtained from MCMC modelling for 13 dark remnant candidate events. Multiple 
entries per event indicate different degenerate solutions found. The full table with the parameters for all 59 events is available online. 
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allax. The mass and distance of a leasing object can be then 
derived from: 



KTTe 


helio^helio 

Mrel 

AvTTe 


Dl 


1 

^helio^helio^j^ + l/Ds 


( 2 ) 

(3) 


where we used the fact that the angular size of the 
Einstein radius can be rewritten as product of the length 
of the vector of the heliocentric relative proper motion 
Iftreil = \^J^L — Msl between lens (L) and source (S) and 
the events time-scale Proper motions and the time-scale 
should be measured in the same frame, either geo- or he¬ 
liocentric. We therefore converted the time-scales of events 
obtained in our geocentric parallax model into the heliocen¬ 
tric frame, following Skowron et al. (2011) (in other words, 
the motion of the Earth was taken out from the relative 
proper motion computation, hence if the lens was observed 
from the Sun, it would have the time-scale of tg®*”, hereafter 
denoted Ie). 


4.1 Proper motions 

Since we do not know the actual relative proper motions for 
our events, we had to assume it. In order to constrain the 
value of /irei, we assumed the sources belonged to the bulge 
population and had bulge-like distribution of velocities, and 
lenses are in the disk and follow disk’s velocity and density 
distributions. We could then draw a random relative motion, 
composed of random source and lens motions, and then use 
the prior on how likely is that random relative proper mo¬ 
tion given standard velocity and stellar density distributions 
for the bulge and disk (e.g.,Batista et al. 2011). The val¬ 
ues were slightly adjusted to match the recently measured 
proper motions by Calamida et al. (2014) based on HST (see 
below). 

Based on the CMD positions (Fig. 2) of the sources we 
assumed that sources belong to the bulge population and 
are located at a fixed distance of 8 kpc. This assumption is 
also in agreement with the results of the simulation of NS 
and BH microlensing events done by Oslowski et al. (2008). 
The lenses, on the other hand, are distributed at various 
distances within the Galactic disk. The closer the lens, the 
more likely a measurable parallax effect will be produced, 
hence we expect a bias towards nearby lenses. Therefore, we 
assumed here bulge-like motions for the sources and disk-like 
motions for the lenses. Because the direction of the vector 
tte coincides with the direction of the relative proper motion 

we were able to narrow the selection of valid combina¬ 
tions of ps and fiL using the angle defined by tte. 


4.2 Probability Density Functions 

MCMC parallax models for each event returned samples of 
solutions (sets of all microlensing parameters with paral¬ 
lax), whose densities reflect their probability. Mixing those 
values with lens and source proper motions, we obtained 
probability densities for masses and distances of our events, 
using Equations (2) and (3). Note that any parallax model 
very often had two or more equally likely separate solutions 
(e.g., Gould 2004, Smith et al. 2005), therefore we carefully 


investigated MCMC results for each event to make sure the 
entire parameter space was explored and all possible solu¬ 
tions were discovered. 

In order to transform the MCMC samples in microlens¬ 
ing parameter space to probability density functions (PDFs) 
for mass and distance, we transformed their density using 
the Jacobian and priors on the relevant parameters, follow¬ 
ing Batista et al. (2011). In our case, however, the relative 
proper motion (/irei) was one of the independent variables 
on the side of the microlensing parameters (as we provided 
its value in order to compute the mass). Each MCMC sam¬ 
ple was weighted with the following weight on microlensing 
rate (F): 


dtEdteM^VTE = :i^K*,J/,^)/(l‘reO[fl(M)M] 


Dl0e\ Ie 
te / tte 


where v is the stellar density distribution, /(p7ei) is the 
distribution of relative proper motions, g{M) is the mass 
function. For the stellar distributions we assumed a thick 
Galactic disk with 0.6 kpc and 2.75 scale height and length, 
respectively. Velocities of the bulge sources in the galactic 
coordinates were assumed as (0,0) km/s in the heliocentric 
frame with a dispersion of (80,80) km/s. For the disk we 
assumed (0,200) km/s and a dispersion of (40,55) km/s. The 
mass function was assumed to be oc 

The individual PDFs for masses of all 13 dark rem¬ 
nant candidates are shown in Fig. 4. Only shown are mass 
distributions for the solution with the highest dark remnant 
probability. Lower panel in Fig. 4 shows histogram of median 
masses of remnant candidates (solid line). The median mass 
distribution corrected for the detection efficiency is drawn 
with dashed fine. Table 5 shows medians and 1 sigma errors 
(for masses and distances found for all Red Clump source 
lenses, i.e., the events for which those parameters were pos¬ 
sible to estimate. The median values were computed from 
1-D likelihood distributions of each parameter marginalized 
over all other dimensions. The dark remnant candidates are 
included in that table. 


4.3 Dark Remnants Candidates 

Stellar remnants acting as lenses should be dark, i.e., they 
should not contribute to the overall observed light. Mi¬ 
crolensing allows us not only to constrain the mass of the 
lens, but also its light output. 

Having drawn a relative proper motion for each sample 
in the MCMC solution, we computed the equivalent mass 
and distance of the lens, following Eqs. (2) and (3). Then, 
for a given mass and distance, the luminosity of the lens was 
derived, using a luminosity-mass relation for main sequence 
stars: L ~ M'^ for masses less than 2 Mq and L ~ 1.5 x 
for M > 2Mq. 

The computed brightness of the lens was then compared 
with the amount of additional light (blended light) in the 
microlensing event, after subtracting the light of the lensed 
source, which is returned in the microlensing model as fs- 
Note that in the crowded Bulge region, blending occurs even 
without microlensing, when the stars are closer to each other 
than a typical seeing radius (~1 arc sec in case of OGLE ref¬ 
erence images), meaning that in case of microlensing events 
there will often be yet another (third) source of light, not 
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Mass [Mq] 


Figure 4. Probability density functions (PDF) for masses of 13 
dark remnant candidates (red thick lines) as obtained for parallax 
microlensing events. Only the highest dark remnant probability 
solutions are shown. Thin vertical blue line for each PDF indicate 
its median value. Bottom panel shows a normalised histogram of 
median masses as measured (solid line) and corrected for detec¬ 
tion efficiency (dashed line). 



related to the microlensing event. Hence, by attributing the 
entire blending light to the light of the lens, we place an 
upper limit on the flux originating from the lensing object. 
The brightness of the blend was corrected for the interstellar 
extinction using parameters from Nataf et al. (2013). The 
bottom panel in Fig. 5 shows a comparison between the ex¬ 
pected brightness of the lens given its mass and distance and 
modelled blending light for a Main Sequence lens, OGLE3- 
ULENS-PAR-42. 

We define the probability that the lens is a dark rem¬ 
nant as the ratio of the integrated density when there is not 
enough blending light to account for a luminous Main Se¬ 
quence lens with a given mass to the total density, i.e.,the 
integrated area in the bottom-left triangle in Fig. 5. 

Table 4 lists all parallax events studied here and in¬ 
cludes the probability for each of the events being a dark 
remnant for each solution. Note, some events had more than 
2 degenerate solutions and probabilities were derived for 
each of them separately, as they yielded different blending 
parameters, masses and distances. 


5 DISCUSSION 

There has been more than 15,000 microlensing events de¬ 
tected so far since the beginning of the microlensing surveys 
in early 1990s. Wyrzykowski et al. (2015) compiled a list 
of 3,500 pure standard microlensing events from OGLE-III 
data (2001-2009), the largest homogeneous sample of stan¬ 
dard events. Those events prove useful for statistical studies 
of mass distribution in the Galaxy, however, they are not 
suitable for discovering individual remnant lenses since they 
provide only one parameter used in mass determination (ts). 


Figure 5. OGLE3-ULENS-PAR-42, an example of a non¬ 
remnant lens parallax event. Top: mass-distance probability den¬ 
sity. Bottom: comparison of computed blended light and lens light 
expected for a given mass and distance. Shown is the solution 
with tto < 0 and probability 25.6 per cent for dark remnant (inte¬ 
gral below dashed line). The lens is most likely a ^ 0.4A4q Main 
Sequence star at Ri 2 kpc. 

Here we have carefully selected those long events which ex¬ 
hibited parallax signatures due to Earth motion, and were 
free of any other additional effects, (e.g., binary lens), as ro¬ 
bust measurements of tte guarantees a good mass estimate. 

By studying the light curves of 150 million sources mon¬ 
itored frequently for 8 years by the OGLE-III project we 
have found 13 candidates for dark stellar remnants acting 
as lenses in microlensing events, with probabilities higher 
than 75 per cent. There would be 19 candidates given a 50 
per cent probability threshold and 5 for a conservative 90 
per cent threshold. 

For each candidate the mass and distance of the lens 
plus its light, to exclude regular stellar lenses, were esti¬ 
mated based on the microlensing parallax model fit to the 
light curve (for parameters Ie and tte), and the unknown rel¬ 
ative proper motion was drawn from observed empirical dis¬ 
tributions. We therefore obtained probability density func¬ 
tions for the mass and distance of each remnant, as well as an 
accumulated probability that a given lens is non-luminous. 

5.1 Proper Motions and Distances 

The relative proper motion (PM, or ^rei) between the source 
and the lens is the key missing ingredient for an exact de- 
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Table 4. Red Clump source events and their probabilities of hav¬ 
ing a dark remnant lens. Probabilities are computed for all solu¬ 
tions and events are sorted according to the highest probability 
of any of the solutions. PAR-08 is the last one included in the 75 
per cent threshold. 


Lens name 

OGLE3-ULENS- 

UQ > 0 
■^EN > 0 

tiQ > 0 
Tl’EN < 0 

tiQ < 0 
■^EN > 0 

tiQ < 0 
TTEN < 

PAR-02 

99.2 

_ 

_ 

99.8 

PAR-13 

- 

99 

97.9 

- 

PAR-05 

95.4 

- 

- 

94.1 

PAR-07 

- 

92.1 

85.9 

- 

PAR-04 

91.2 

- 

- 

91.7 

PAR-19 

- 

89.2 

- 

62.6 

PAR-33 

87.9 

41.6 

22 

60.6 

PAR-15 

- 

87.1 

73.9 

- 

PAR-09 

- 

86.9 

- 

- 

PAR-39 

- 

86.8 

- 

85.9 

PAR-28 

77.8 

- 

- 

85.7 

PAR-27 

85.3 

74 

52.8 

60 

PAR-08 

76.6 

- 

- 

72.5 

PAR-24 

71.6 

45.3 

51.4 

41.8 

PAR-34 

63.6 

- 

70.8 

- 

PAR-03 

65.2 

- 

- 

61.5 

PAR-22 

57.2 

- 

- 

- 

PAR-01 

- 

- 

53.8 

- 

PAR-44 

- 

53.2 

50.3 

- 

PAR-06 

35.7 

- 

- 

43.6 

PAR-42 

- 

43.2 

25.6 

- 

PAR-36 

- 

42.5 

- 

32.1 

PAR-16 

- 

33.4 

40.9 

- 

PAR-23 

- 

38.7 

27.8 

- 

PAR-47 

- 

33.3 

- 

- 

PAR-14 

5.5 

- 

- 

12.0 



termination of the mass and distance of the lens. Therefore, 
we have used the prior on distribution of PM measured for 
Galactic bulge and disk stars (Calamida et al. 2013), as¬ 
suming our sources are in the bulge and the lenses in the 
disk. 

However, if a source and lens were both located in the 
Galactic disk, they might together have very small value of 
Prei if their projected motion is almost parallel. Such PM 
would decrease the measured mass of the lens (Eq. (2)) and 
would place the lens at a further distance (Eq. (3)), hence 
increasing a chance the lens is a normal disk star. 

Since the microlensing phenomenon amplifies the light 
of a distant source star, it is possible to determine the true 
brightness of the source, despite severe blending conditions 
typical for the bulge region (see e.g.,Wozniak & Paczyhski 
1997, Wyrzykowski et al. 2006, Smith et al. 2007). Our sam¬ 
ple of parallax events (59) was limited to only those with a 
source located in the Red Clump Region (26 events). Since 
the vast majority of Red Clump Giants towards the bulge are 
located at about 8 kpc, we assumed that our RCR sources 
have Ds = 8 kpc and they follow the PM distribution of the 
Bulge stars. There are of course rare cases of RC stars in 
the Galactic disk, but they would be shifted on the CMD as 
brighter and less affected by extinction, hence might fall out 
of our RCR box. The same applies (but with the inverse ar¬ 
gument about CMD location) to even rarer examples of Red 
Clump stars from the far side of the Galactic bulge. Note 
that our mass determination (Eq. (2)) does not depend on 
the assumption on Ds, but only on firei- The distance to the 
lens, computed from Eq. (3), for a typical time-scale in our 
sample of ~200 days and parallax tte ~ 0.05, can vary by 
at most 0.5 kpc if Ds is different from 8 kpc by ±1 kpc. 


Table 5. Most probable (1-D median and 1 sigma) masses, dis¬ 
tances and total blend luminosities for all events with Red Clump 
sources. Multiple entries are shown for different solutions found, 
indicated in the last column. 
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for small values of fi^ei ~ 2. The error-bar from the MCMC 
modelling for Dl is similarly wide. 

Another crucial assumption is that the lenses are lo¬ 
cated in the disk and follow the PMs of the disk stars. How¬ 
ever, the events we use they exhibit a clear microlensing 
parallax signal and a large parallax signal indicates a nearby 
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lens. For our assumed Ds = 8 kpc and typical ttb ~ 0.05 
most lenses should be located at Dl < 6 kpc, hence clearly 
lying in the Galactic disk. On the other hand, a lens located 
in the Bulge would have a small and hard to detect tte. A 
relative PM distribution for such microlensing event would 
not differ much from the bulge-disk one, but even an in¬ 
creased chance for a smaller /Trei would be compensated by 
small value of tte and the mass measurement would not be 
affected significantly. 

Of course, our remnant lens candidates, especially black 
holes, do not have to follow general disk motions at all - for 
example they might reside in the Galactic halo or be mov¬ 
ing as a result of a significant birth kick (Fryer et al. 2012). 
However, in such situations, their relative motion with re¬ 
spect to the source star would be even larger than assumed 
for the disk. Hence, our disk motion assumption provides a 
safe lower limit on the mass of the remnant. 

We have to note that the PDFs for the mass mea¬ 
surements are still wide and range often over an order of 
magnitude. If the assumptions on relative proper motions 
are invalid in case of some events, their masses could dif¬ 
fer significantly. In principle, the lower mass end is heavily 
populated by white dwarfs, the distribution of which peaks 
at around 0.6 Mq (Kepler et al. 2007). Most white dwarf 
lenses with absolute magnitude in range between 2-15 mag 
could still fit within the light of the blend of most of our 
candidates, at least in the tails of the probability distribu¬ 
tions. We note that our mass distribution does not show any 
prominent peak at masses where a vast majority of white 
dwarfs should be. Our experiment is less sensitive for lower 
masses (see Fig. 3) and the detection efficiency of our search 
pipeline (which does not include expected Galactic remnants 
population densities) reflects that (see dashed line in Fig. 4, 
however, there still should be significantly more WDs than 
neutron stars and black holes. This indicates that there still 
could be a strong bias in parallax microlensing events for 
heavier lenses and our search is still not completely sensi¬ 
tive to white dwarf lenses. 


5.2 Individual events 

Below we discuss some of the individual microlensing events 
having a candidate dark remnant lens. 

OGLES- ULENS-PAR-01 

OGLE3-ULENS-PAR-01 is the event with the longest dura¬ 
tion and is probably the longest event ever observed lasting 
2840.8 days (~7.8 years). The event was not found by the 
real-time OGLE-HI data analysis (EWS), probably because 
its deviation had started already in 2002. In its light curve 
(shown in Fig. 6) it exhibits multiple peaks due to the par¬ 
allax effect modulation, making PAR-01 the second multi¬ 
peak microlensing event after OGLE-1999-BUL-19 (Smith 
et al. 2002). The large amplitude parallax signal immedi¬ 
ately indicates that the lens is nearby. There was only one 
physically correct solution for parallax (uo < 0) returned 
by the MGMC modelling. MCMG modelling corrected for 
Galactic priors returned the most likely distance of just 1.3 
kpc and mass of 1.0 Mq. At such a mass and distance and 
with a blending parameter of fs = 0.68 there was some- 
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Figure 6. Light curve and standard (green) and parallax (red) 
microlensing models and their residuals for the longest lasting 
multi-peak event OGLE3-ULENS-PAR-01. There was only one 
solution found and it has an equivocal 53.8 per cent probability 
for a 1.0 Mq dark lens at about 1.3 kpc. The width of the parallax 
model curve (red) indicates 1 cr uncertainty in the model. 

what a low probability (53.8 per cent) that the lens is non- 
luminous. At this moment we can not rule out between a 
main sequence star or a neutron star, however, in the tail 
of the probability distribution there is still a possibility of 
a massive white dwarf, if the relative proper motion be¬ 
tween the bulge source and the nearby lens was significantly 
smaller than few mas/yr. On the other hand, since the lens 
is near (consistent with the large parallax signal), its proper 
motion could be different than that assumed for standard 
disk stars. If significantly higher, then the mass of the lens 
would also become significantly heavier, hence it would be 
more difficult to explain the blended light with a nearby 
massive main sequence star and a nearby neutron star and 
even a black hole solution could be considered. 

Further observations of this object, including high reso¬ 
lution imaging, could potentially help resolve the source and 
the lens and a non-detection of a lens could indicate another 
dark remnant. 

OGLES- ULENS-PAR-02 

OGLE3-ULENS-PAR-02 is the second longest event in our 
sample, with a length of 2306.4 days (6.3 years), spanning 
until the very end of the OGLE-HI phase in 2009. A small 
degree of amplification is also still visible in the OGLE-IV 
data from 2010. The light curve with the most probable 
microlensing model of PAR-02 is shown in Fig. 7. 

The MCMC modelling yielded two degenerate solutions 
for Uo > 0 and uo < 0 (see Table 3), with heliocentric time- 
scales of tE = 296^.7 days and Ie. = ShOlg days. Parallax 
parameter tte solutions are shown in Fig. 8. 

The light curve for this event is shown in Fig. 7 along 
with the standard and parallax microlensing models. The 
deviation from the standard model, shown at the bottom, 
has an amplitude of only ±3 per cent, however, thanks 
to frequent monitoring and good photometry, the parallax 
anomaly is clearly visible with yearly fluctuations and is well 
reproduced by the model. 
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Table 6. Parallax model parameters of OGLE3-ULENS-PAR-01, 
the longest parallax microlensing event with multiple peaks in its 
light curve. 
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Figure 7. Light curve and standard (green) and parallax (red) 
microlensing models and their residuals for the event OGLE3- 
ULENS-PAR-02. The solution shown (ug < 0) has a 99.8 per 
cent probability of a dark lens of 8.7Mq lens at 1.8 kpc. 
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Figure 8. Parallax degeneracy found in the photometric data 
parallax models for OGLE3-ULENS-PAR-02 event. Blue, green 
and red contours indicate 1,2,3 sigmas, respectively. 


The two degenerate solutions differed significantly in 
the computed relative source light contribution to the over¬ 
all light: /s = l.Ollg gg and /s = 0.64 ± 0.04, with a slight 
preference for the fs = 1.01 solution (that preference cor¬ 
responding to a difference in of = 8.5). A blending 
parameter close to 1.0 means that the source dominates the 
light at the place of the event and indicates very little room 
for any extra light from the blend, in particular from the 
lens. The distribution of allowed blend brightness in our so¬ 
lution is shown in the lower panel of Fig. 9. The blend can 
still have a brightness of about 17.7 mag, but for our most 
likely mass and distance, this most likely is not the lens 
itself. 

OGLE3-ULENS-PAR-02 has the highest probability for 
being a dark remnant among our candidates - both paral¬ 
lax solutions indicate probabilities above 99 per cent that 
the lens is a dark remnant. The somewhat more probable 
of the two solutions (uo < 0) gives the lens at a distance 
of 1.8 kpc and with a mass of 8.711 '^Mq, but the other so¬ 
lution (with some blending) yielded an even higher mass of 
9 .3l4'3M0 at 2.4 kpc. PAR-02 is the heaviest of our dark 
remnant candidates, hence it is most likely an isolated black 
hole. Interestingly, this black hole candidate is the only one 
in our sample where the lens mass is above the mass cut-off 
observed in black holes in X-ray binaries (Ozel et al. 2010). 
On the other hand, the detection efficiency in our pipeline 
for microlensing events with parameters similar to PAR-02 
is relatively low (see Fig. 3), since the smaller value of tte, 
the more massive the lens. Hence there might be more black 
hole lenses with masses above 6 Mq, which do not exhibit 
strong parallax signals in our data and are hidden among the 
standard microlensing events. Fig. 9 shows the probability 
density for mass and distance for this solution as well as 
a comparison between computed blend light with expected 
luminosity of the lens for its mass and distance. 

The mass-distance-blend light derived for the other so¬ 
lution (uo > 0 and/s = 0.64 ± 0.04) also indicates a dark 
remnant/black hole with very high probability (99.2 per 
cent). There is blended light in this solution, however, it 
is most likely to be coming from an unrelated nearby blend¬ 
ing star since a main sequence star with mass of 9.3 Mg at 
distance of 2.4 kpc should be significantly brighter. The to¬ 
tal blend light is estimated at /giend ~ 16.3 mag (extinction 
corrected). 

If, nevertheless, the lens was indeed the blend in this 
solution, given its brightness and assuming again a Main 
Sequence star from the disk, the lens would need to have a 
mass in the range between 0.9 and 1.1 Mg and then the rela¬ 
tive proper motion would need to be about 0.9 mas/yr. Such 
a small value of proper motion would mean that both source 
and lens move almost in parallel (projected). This is not im¬ 
possible, but the chance is extremely small for bulge and disk 
stars. Such an interpretation can be tested with high resolu¬ 
tion Adaptive Optics (AO) observations of the system after 
10 or more years, when the source and lens would have sepa¬ 
rated by «10 mas (e.g., Kozlowski et al. 2007, Pietrukowicz 
et al. 2012, Batista et al. 2015). With both source and lens 
being relatively bright, a deformation in the point-spread 
function should allow us to measure such delicate blend. On 
the other hand, a non-detection of lens light would allow us 
to rule out this solution with high confidence. 

Since both solutions indicate a massive lens, the event 
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Figure 9. OGLE3-ULENS-PAR-02, the most massive and most 
probable isolated black-hole candidate. Top: mass-distance prob¬ 
ability density. Bottom: comparison of computed blended light 
and lens light expected for a given mass and distance. Shown is 
solution with uq < 0 and probability 99.8 per cent for dark rem¬ 
nant (integral below dashed line). The lens light axis stretches to 
the very bright magnitudes. 

OGLE3-ULENS-PAR-02 is the best black hole candidate in 
our sample, with its mass of 8.7 or 9.3Mq, depending on 
solution, and strong constraints on the lens light. Neverthe¬ 
less, further detailed AO and X-ray follow-up is necessary 
in order to rule out other solutions and to potentially detect 
the BHs interaction with the interstellar medium (e.g., Agol 
& Kamionkowski 2002, Fender et al. 2013). 

OGLES- ULENS-PAR-05 

OGLE3-ULENS-PAR-05 is another very strong candidate 
for a dark remnant lens. Its lightcurve, shown in Fig. 10, 
shows a very long-term and strong parallax; anomaly, allow¬ 
ing for a good measurement of ttb. There were two solutions 
found. Both solutions have blending /s very close to 1, in¬ 
dicating no blending and leaving no room for light from a 
massive lens. Probabilities for dark remnant lens were above 
95 per cent for both solutions. The most probable mass of 
the lens was estimated at 3.3_]^'5 and 4 :. 8_2 ^M 0 , for positive 
and negative uq solutions, respectively. Figure 11 shows the 
mass-distance probability density and comparison of com¬ 
puted blended light and lens light. 

If the mass of this dark remnant was indeed as our most 
likely estimate, this object would be residing in the very 



Figure 10. Light curve and standard (green) and parallax (red) 
microlensing models and their residuals for the event OGLE3- 
ULENS-PAR-05. The solution shown (uq > 0) has a 95.4 per 
cent probability of a 3.3 Mq lens at 2.9 kpc. 

gap between known neutron stars and black holes. Kiziltan 
et al. (2013) suggested that the mass distribution of NSs, 
despite its apparent concentration between 1 and 1.5 Mq, 
can also spread to higher masses, also above 2 Mq. On the 
other hand, the mass distribution for black holes, found in 
X-ray binaries, have a cut off at about 6 Mq (Ozel et al. 
2010). Microlensing parameters of OGLE3-ULENS-PAR-05 
place it in the region of relatively low detection efficiency 
(see Fig. 3). A detection of OGLE3-ULENS-PAR-05 indi¬ 
cates there might exist a population of isolated NSs or BHs 
present in the mass distribution gap, in agreement with cur¬ 
rent theoretical predictions (Fryer & Kalogera 2001). 

5.3 Mass distribution 

Fig. 4 shows individual mass probability density functions 
for the most probable remnant solution of our 13 events with 
remnant lens candidate. Median masses, listed also in Table 
5, are marked with blue vertical line and their distribution 
is shown in the lower panel of Fig. 4. It can be seen that 
correcting for the detection efficiency only mildly modifies 
the shape of the histogram, with the lowest and the high¬ 
est masses being detected with somewhat lower detection 
efficiency. If indeed all our 13 lenses are due to remnants 
and we use their most probable masses, the derived approx¬ 
imate mass distribution has a peak between 1 and 2 Mq, 
most likely due to a mixture of single neutron stars and 
massive white dwarfs. The low-mass end is then populated 
by white dwarfs, however, as noted above, we seem to be not 
complete for WDs. In the high-mass end of the distribution 
we see a continuum of masses, from 1 to 10 Mq. Strikingly, 
the shape of this mass function resembles the synthetic one 
obtained by Fryer & Kalogera (2001), with a steady decline 
of masses of neutron stars and black holes, with no mass gap 
between the two populations of compact remnants. 

5.4 Astrometric microlensing 

Microlensing events with heavy and relatively nearby lenses 
will have large angular Einstein radii, 0e, where 0e oc 
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Figure 11. OGLE3-ULENS-PAR-05, a dark remnant lens candi¬ 
date from the mass gap. Top: mass-distance probability density. 
Bottom: comparison of computed blended light and lens light ex¬ 
pected for a given mass and distance. Shown is solution with 
uq > 0 and probability 95.4 per cent for dark remnant (integral 
below dashed line). 

and Ttrei = AU{l/Di^ — l/Ds). As shown in, 
e.g., Walker (1995), H0g (1995) and Dominik & Sahu (2000), 
the motion of unresolved images in a microlensing event pro¬ 
duces a characteristic displacement of the observed centroid 
of light of the source star of an order of 9^- An astromet¬ 
ric deviation from the unperturbed motion of the source is 
computed as S{u) = where u is the vector between 

source and lens positions. It reaches its maximum amplitude 
of 5max ~ 0.7 6 e for Wo = \/2. It means that by knowing u 
from photometry and measuring the centroid displacement 
one can immediately derive the value of 6 e for microlens¬ 
ing events without the need of any other additional effects 
(finite source or caustic crossing). Following Eq. (1), having 
measured both tte and 0 e one can compute the mass of the 
lens directly. 

Fig. 12 shows distributions of probable ranges of 9e for 
dark lens candidates. They range from a fraction of mas 
to couple of milliarcseconds. According to de Bruijne et al. 
(2014) ESA’s Gaia mission at the end of the mission should 
achieve astrometric precision better than 1 mas for stars as 
faint as E ~ 19 mag. Gaia is collecting observations since 
mid-2014 and is expected to do so until 2019/2020. There¬ 
fore, for recently on-going microlensing events with heavy 
remnant lenses, similar to our sample, we will have superb 
astrometric time-series provided by Gaia. Once combined 



Figure 12. Distribution of probability for angular Einstein radius 
(0E ) for 13 parallax microlensing events with dark remnant lens. 
Vertical thin line indicate median, most probable value. Most 
events with 6 e > 1 rnas will have their Einstein radii measured 
from Gaia astrometry. 

with ground-based (or space-based) parallax measurements, 
a mass and distance of the lens will be measured without 
the need of any assumptions about proper motions. 

Several searches for those minute astrometric centroid 
displacements are being conducted, however, so far none 
were reported positive. In the near future, however, the 
Gaia space mission (e.g.,de Bruijne et al. (2014)) will pro¬ 
vide sub-milliarcsecond positional time-series for a billion 
stars in our Galaxy, down to E ~ 20 mag. Such precise as¬ 
trometry will be most suitable for detecting and measuring 
microlensing astrometric signals (Belokurov & Evans 2002, 
Wyrzykowski & Hodgkin 2012). Fig. 13 shows examples of 
simulated astrometric trajectories for our best BH candi¬ 
date, along with simulated Gaia observations, as if Gaia was 
launched a decade earlier. Similar events with black hole 
lenses will produce similar paths, which will allow for inde¬ 
pendent measurements of 6 e- When combined with tte (ob¬ 
tained either from Earth orbit accelerations or from space, 
e.g.,from Spitzer). Gaia will provide direct measurements 
of masses of lenses for hundreds of microlensing events. In 
particular, dark lenses would produce easy-to-measure as¬ 
trometric displacements, as the astrometry would be less 
affected by blending. Detections of neutron stars and black 
holes via astrometric and photometric channels will open a 
new path of studying the mass function of dark remnants. 
However, we stress that it would not be possible without 
long-term photometric monitoring programmes. 


6 CONCLUSIONS 

We have found 13 candidates for stellar dark remnants us¬ 
ing microlensing events with measurable parallax effects. 
Among them are both neutron stars and black holes, as 
microlensing can not distinguish between the two, but we 
used the modelled light to rule out stellar luminous lenses. 
We estimated the masses and distances of those objects as¬ 
suming standard proper motions for the lenses. The most 
massive lens found in our sample is OGLE3-ULENS-PAR- 
02 with S.TlI'yM© for the more likely of the two solutions. 
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Figure 13. Simulations of example astrometric paths (solid lines) for a black-hole event similar to OGLE3-ULENS-PAR-02 for pj-el = 4 
mas/yr, but with different directions of that vector. Dashed lines show unperturbed source motion with parallax. Dots show mock Gaia 
measurements with Gaia sampling and expected accuracy. Note that error-bars along the scanning direction (blue) are significantly better 
than across the scanning (thin red lines). 


and 9 . 3 I 4 jMq for the other solution. Another candidate, 
OGLE3-ULENS-PAR-05, with most likely mass of 3.3 or 
4.8 Mq, seems to reside in the mass gap in the mass dis¬ 
tribution for NSs and BHs and indicates there might be a 
population of isolated massive remnants filling the observed 
gap, in line with the current theoretical predictions. 

Our mass estimates rely on certain assumptions on dis¬ 
tances and relative proper motion. In principle, it can not 
be currently fully ruled out that some of our dark remnant 
candidates are simply main sequence stars in the disk, which 
happened to co-move in parallel to a bulge source star, yield¬ 
ing very small relative proper motion. Also, we can not fully 
exclude a possibility that the lower mass lenses (around and 
below 1 Mq) are nearby white dwarfs. However, the most 
heavy lenses found in our study still remain the best can¬ 
didates for isolated black holes within or outside the mass 
gap between neutron stars and black holes. The approxi¬ 
mate mass function obtained for our 13 candidates resembles 
the one obtained in the stellar evolution synthesis and indi¬ 
cates no mass gap above 2 Mq. The sample of microlensing 
events with a robust parallax measurement will be further 
expanded with many events being currently found in the 
OGLE-IV survey. 

We also note that the events with the greatest length 
are more likely to be due to massive lenses and exhibit strong 
parallax effects, allowing us to constrain the mass of the lens. 
Gonducting a search on combined OGLE Bulge data from 
1997 (OGLE-II, OGLE-III and OGLE-IV) would yield dis¬ 
coveries of the most massive lenses and black holes ever. On 
the other hand, the currently on-going OGLE-IV project is 
finding about 2000 events every year. Those events are also 
currently observed by Gaia, which at the end of its five-f 
years long mission will provide sub-milliarcsecond astromet¬ 
ric time-series for them. Combining both ground-based su¬ 
perb photometry and space astrometry will yield many dis¬ 
coveries of black hole lenses among microlensing events and 
will allow us to put constrains on the remnant mass function. 
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